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OO ■ Abstract 

(N : 

A d-dimensional spacetime is "axisymmetric" if it possesses an SO{d—2) isometry group whose 
' orbits are (d— 3)-spheres. In this paper, algebraically special, axisymmetric solutions of the higher 

dimensional vacuum Einstein equation (with cosmological constant) are investigated. Necessary 
^ • and sufficient conditions for static axisymmetric solutions to belong to different algebraic classes 

are presented. Then general (possibly time-dependent) axisymmetric solutions are discussed. All 
axisymmetric solutions of algebraic types II, D, III and N are obtained. 

00 ■ 1 Introduction 

^ ' 1.1 Background 
'NT I 

■r^ij- ■ A d-dimensional spacetime is "axisymmetric" if it possesses an SO{d — 2) isometry group whose 

i orbits are {d — 3)-spheres. There are several motivations for studying axisymmetric solutions of the 

I higher-dimensional vacuum Einstein equation (with cosmological constant) 

> 

X 



Rau = ^Qav (1) 



^ \ These include the problem of finding an exact solution describing a black hole bound to a 3+1 
dimensional brane in the (single brane) Randall-Sundrum model ^j, and determining the phase 
structure of General Relativity with a compactified dimension [2]. 

In d = 4 dimensions, all static axisymmetric solutions of the vacuum Einstein equation (with 
A = 0) were obtained by Weyl, who showed that they are characterized by a single axisymmetric 
harmonic function in E? (see e.g. [3]). Weyl's result has been generalized to higher dimensions: 
the class of solutions of the d-dimensional vacuum Einstein equation (with A = 0) admitting d — 2 
commuting, orthogonal, non-null Killing fields is specified by d — 3 axisymmetric harmonic functions 
in E? |4J. If one of the Killing fields generates time translations and the others generate rotations 
then these solutions have isometry group R x U{lY~^ , generalizing the R x U{1) symmetry of Weyl's 
solutions. However, these solutions are not axisymmetric for d > 4. 

It is desirable to know the general static axisymmetric solution va. d> A dimensions but, unfortu- 
nately, the Einstein equation cannot be solved analytically for (i > 4 (or even for d = A with A 7^ 0). 
The impediment arises from the curvature of S'^~^ [5j. Note that S"^~^ is fiat if d = 4, which is why 
the Einstein equation can be solved for d = A. 

The goal of this paper is to determine whether the Einstein equation can be solved analytically 
if one makes the additional assumption that the spacetime is algebraically special. In d = 4, an 
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algebraically special static, axisymmetric spacetime must be type D (or O). For A = 0, the only 
such solutions are, in Levi-Citiva's evocative terminology, the A-metrics, the B-metrics, and the 
C-metric [6]. The A-metrics are labelled by the parameters k G {—1,0,1} and M ^ 0. The metric 
takes the generalized Schwarzschild form 

ds^ = -Uir)dt^ + ^^+rW„ (2) 

where U{r) = k — 2M/r, and dil"^ is the metric on a space of constant curvature with sign k. The 
B-metrics are analytic continuations of the A-metrics in which the time coordinate t is Wick rotated 
to a spatial coordinate (j) and dClf. to a Lorentzian metric of constant curvature 

ds^ = U{r)d^'^ + ^ + r^dS^. (3) 

The C-metric describes a pair of black holes being accelerated apart by a conical singularity [7]- 

The first method for algebraic classification of the Weyl tensor in higher dimensions was proposed 
by de Smet [8J. His spinorial method works only for the case d = 5. We shall follow the alternative 
algebraic classification of Coley, Milson, Pravda and Pravdova (CMPP) [9], which applies for general 
d. This method uses "Weyl aligned null directions" (WANDs), which generalize the concept of 
principal null directions (PNDs) used for algebraic classification in d = 4 dimensions. We shall 
explain the classification scheme in more detail below. For now, we just state that a > 4 spacetime 
need not admit a WAND, in which case it is algebraically general (type G), and that there is a notion 
of a "multiple WAND" that generalizes the 4d notion of a repeated PND. If a WAND exists, but 
not a multiple WAND, then the spacetime is type I. If a multiple WAND exists then the spacetime 
is type II, D, III, N or O. The conditions for these more special algebraic types will be described 
below. Finally, we note that WANDs need not be discrete for d > 4, and the terminology PND is 
reserved for the case in which there are finitely many WANDs. 

An important result in 4d is the Goldberg-Sachs theorem [6]. This states that, for a solution of 
the Einstein equation ([1]) that is not type O, a null vector field is a repeated PND if, and only if, it 
is tangent to a shear-free null geodesic congruence. This is not true in more than four dimensions 
[101 lllj : a multiple WAND may not be geodesic, or it may be geodesic but shearing. An example 
of the former behaviour (with A > 0) is the direct product dS^ x S'^^^ (for d > 4), which is type 
D [12]. Any null vector field in dS^ is a multiple WAND. Obviously not all such vector fields are 
geodesic. An example of the latter behaviour is the black string solution given by the product of 
the 4d Schwarzschild solution with a fiat direction, which is also type D yL2j. The multiple WANDs 
are the repeated PNDs of the Schwarzschild solution. These are geodesic (by the Goldberg-Sachs 
theorem). The associated congruences of geodesies expand in the Schwarzschild directions but not 
in the flat direction, hence they are shearing0 

A partial generalization of the Goldberg-Sachs theorem does hold for vacuum spacetimes of type 
III or type N, for which it has been shown that the multiple WAND must be geodesic [lljl^ 



1.2 Summary of results 

In this paper, we start (in section [2|) by considering static, axisymmetric solutions and determine the 
condition for them to be algebraically special. It was shown in Ref. [12] that a static solution must 
be of algebraic type G, I, D or O. We derive simple necessary and sufficient conditions for a solution 
to belong to the various algebraic types. We also show that many analytic solutions are type G in 
one open subset of the spacetime and type I in another. This suggests that distinguishing between 

^ See section 5.4 of Ref. [13] for discussion of constraints on the shear of a geodesic multiple WAND. 
^ The argument of Ref. [11] assumes A = but it is straightforward to generalize to A 7^ 0. 



2 



type G and type I solutions is not very useful in practice, and that the type I condition alone will 
not be much help in finding new solutions (just as in 4d, where type I is algebraically general) H 

In the rest of the paper, we relax the condition of staticity and consider general (possibly time- 
dependent) algebraically special axisymmetric solutions. The starting point of our analysis is the 
observation that, for d > 4, the action of S0{d — 2) on S'^~^ must be orthogonally transitive [6], i.e., 
spacetime is locally a product M3 x 5*^"^ with warped product metric 

ds'^ = gab{x)dx^dx^ + E{xfdQ'^, (4) 

for some 3-metric gab and function E(x) on M3, where is the metric on S'^~^ normalized to unit 
radius. The analysis naturally divides into two cases depending on whether or not the WAND is 
axisymmetric, i.e., invariant under SO{d — 2). 

In section [3l we consider the case in which the WAND is axisymmetric. In this case, there 
is a partial generalization of the Goldberg-Sachs theorem. We show that the only axisymmetric 
solution of type II or D with a non-geodesic axisymmetric multiple WAND is dSs x S"^'^ with 
A > 0. For any other axisymmetric type II or D vacuum solution, an axisymmetric multiple WAND 
must be geodesic. As already noted, this is also true for type III and type N pJLj. (Hence any 
axisymmetric vacuum solution with an axisymmetric but non-geodesic multiple WAND must be 
type O, or dS^ x S'^~^.) The axisymmetry implies that the null geodesic congruence tangent to the 
WAND has vanishing rotation, hence it is hypersurface-orthogonal. We determine all solutions with 
an axisymmetric geodesic WAND without assuming a particular algebraic type. The solutions are 
all type II or more specialQ There are several classes. 

• Type O (conformally flat) solutions. Irrespective of axisymmetry, the only such solutions are 
Minkowski, de Sitter, and anti-de Sitter spacetimes. 

• The Schwarzschild solution, generalized to allow for flat or hyperbolic slices (i.e., higher- 
dimensional analogues of the A- metrics) and a cosmological constant. The metric is given 
by equation (|64p . The solution is type D. The null congruence tangent to the WAND has 
vanishing shear and non-vanishing expansion, so these solutions are a subset of the higher- 
dimensional Robinson- Trautman family of solutions (defined to be solutions admitting a null 
geodesic congruence with vanishing shear and rotation and non-vanishing expansion) obtained 
inRef. [IS]. 

• "Black string" solutions obtained, for A = 0, by foliating Minkowski spacetime with (d — 1)- 
dimensional Minkowski, or de Sitter, slices and replacing the slices with a Schwarzschild, 
or Schwarzschild-de Sitter, metric respectively. In the former case, this gives the familiar 
Schwarzschild black string solution. There is an analagous construction for A > based on 
a de Sitter foliation of de Sitter spacetime, and for A < based on Minkowski, de Sitter, or 
anti-de Sitter foliations of anti-de Sitter spacetime. The latter includes the anti-de Sitter black 
string of Ref. [16j. The metric of these solutions is given by equation ()58p . They are all type 
D. The null congruence associated with the WAND has non-vanishing expansion and shear. 

• For A > 0, dSs x S"^'^ is type D. As discussed above, a general multiple WAND of this 
spacetime is non-geodesic. However, any null geodesic congruence in dSs defines a geodesic 
multiple WAND, which is why this solution is mentioned here. Such a congruence may be 
expanding and shearing or non-expanding and non-shearing (in the latter case the solution is 
a special case of the Kundt solutions discussed next). 

^This also suggests that it might be convenient to change nomenclature and reserve the term "algebraically special" 
for spacetimes admitting a multiple WAND, just as in 4d. 

* For A = 0, this implies that these solutions belong to the class of spacetimes discussed in Ref. [14], i.e., those 
admitting a hypersurface orthogonal geodesic multiple WAND. The dependence on the afiine parameter along the 
geodesies was determined in that paper. 
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• Axisymmetric Kundt solutions (section I3.3.4p . A Kundt spacetime is a spacetime admitting 
a null geodesic congruence with vanishing expansion, rotation and shear [6j. Such solutions 
are type II, or more special, for any d > 4 [TTj. In general they involve arbitrary functions of 
time. In our axisymmetric case, these solutions are expressed in terms of solutions of certain 
ODEs that cannot be solved analytically in general. We show that some of these solutions are 
type D or N (but not III). The type N solutions can be obtained analytically. They describe 
gravitational waves in Minkowski (eq. (f99l) ). de Sitter (eq. (jlOSp ) or anti-de Sitter (eq. (11041) ) 
spacetime. The general type D solution is cohomogeneity-1 with surfaces of homogeneity 
M2 X S'^~^ where M2 is 2d Minkowski or (anti-) de Sitter spacetime: 

ds'^ = dz'^ + A{zfdJ:'^ + R{zfdn'^, (5) 

where dS^ is the metric on M2. The functions A{z) and R{z) can be determined analytically 
only in special cases e.g. the product spaces dS^ x S'^^^, dS2 x S'^^'^ and AdS2 x H'^^^, or 
for flat M2 with A = 0. Various solutions of this form have been discussed previously in the 
literature. For positively curved M2, if one analytically continues to Riemannian signature 
(so that M2 becomes 5"^) then these are metrics of the form discussed by Bohm. He proved 
that, for low enough d > 4, and A > 0, there exist infinitely many Einstein metrics on spheres 
and products of spheres of the form ([5]) [18]. for A < he constructed complete, non-compact 
metrics of the form ([5]) pL9j . The Lorentzian interpretation of some of the Bohm solutions has 
been discussed in Ref. [20]. Some singular solutions with flat M2 and A = were discussed in 
Ref. [21], analagous solutions with non-flat M2 were discussed in Ref. [22]. A regular solution 
with d = 5, flat M2 and A < describes the metric dual to the ground state of = 4 super 
Yang-Mills theory on i? x S'"'^ x S''^ (with fermions periodic on S^) [23] . A solution with d = 5, 
negatively curved M2 and A < describes the bulk near-horizon geometry of an extremal 
charged Randall-Sundrum black hole, or the metric dual to the ground state of = 4 SYM in 
AdS2 X 5^ [24j . More generally, solutions of the form (jS]) with A < that are asymptotically 
locally anti-de Sitter presumably describe the metrics dual to the ground state of a CFT in 
M2 X S'^-^. 

For d = 4, it was proposed in Ref. [25] that some type II Kundt solutions describe gravitational 
waves propagating in a "background" spacetime described by a type D Kundt solution. The 
same is true for d > 4: given a type D background of the form ([5]), one can construct explicitly 
axisymmetric type II Kundt solutions which describe gravitational waves propagating along 
the space M2 of this background. (Some particular examples of such solutions were obtained 
in Refs. [26l [2^.) For A < 0, some of these solutions will be asymptotically locally AdS, 
and will describe metrics dual to certain CFT states in M2 x S'^~^ for which there is a null 
energy-momentum flux along M2. 

The second case to consider is when the WAND is not axisymmetric. Acting with SO{d — 2) 
generates a continuously infinite family of WANDs, which suggests that the solution should have an 
enhanced symmetry. This is indeed the case: assuming that the solution admits a multiple WAND, 
we are able to show that SO{d — 2) is enhanced to the de Sitter symmetry 5*0(1, d — 2), with dSd-2 
orbits, and the only non-trivial (i.e. not type O) solutions are: 

• For any A, a Kaluza-Klein bubble solution [28j (i.e. a higher-dimensional analogue of the 
B-metrics) obtained by analytic continuation of the Schwarzschild solution (so that 5"^~^ 
dSci-2)- The metric is given by equation (jl22p . This is type D. 

• For positive A there is a dS(i-2 ^ S'^ solution. This is also type D. 

For both solutions, any null vector field tangential to the dSd-2 orbits of SO{l,d — 2) is a multiple 
WAND. Hence, these are further examples of type D vacuum solutions for which multiple WANDs 
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need not be geodesicJl 

Combining these results, we learn that any algebraically special axisymmetric solution that is not 
encompassed by our analysis must be type I and such that every WAND is either not invariant under 
SO{d—2) or is invariant under S0{d—2) but is not geodesic. In particular, any axisymmetric solution 
admitting a multiple WAND is one of the solutions listed above. It is convenient to summarize our 
results according to algebraic type: 

• Type O: the only type O Einstein solutions are Minkowski or (anti)-de Sitter spacetime. 

• Type N: the only axisymmetric solutions are the type N axisymmetric Kundt solutions. 

• Type III: there are no axisymmetric type III solutions. 

• Type D: all axisymmetric solutions are contained in the following list: Kaluza-Klein bubble; 
dSd-2 X S*^; generalized Schwarzschild; generalized black string; solutions of the form ([5]). 

• Type II: the only axisymmetric solutions are the type II axisymmetric Kundt solutions. 

• Type I: if a WAND of an axisymmetric type I solution is axisymmetric then it is non-geodesic. 

Note that the type D solutions all have isometry groups larger than the SO{d — 2) that was assumed 
initially. 

We can compare these results to those of de Smet, who classified static^ axisymmetric d = 5 
spacetimes belonging to classes 22 and 22 in his classification scheme for A = [8] and A 7^ [29]. 
For d = 5, our list of type D solutions is very similar to the set of the solutions that he foundl§ One 
significant difference is for A = 0, where he found a solution that is not on our list (eq. 4.6 of Ref. 
[8], a "homogeneous wrapped object"). The results of section [2] below show that this solution is type 
G in the CMPP classification. Curiously, no analagous solution with A 7^ was obtained in Ref. 
[29]. Some of the generalized "black string" solutions that we found (eq. ()58p ) do not appear in de 
Smet's results. It would be interesting to understand the connections between the de Smet scheme 
and the CMPP scheme. 

It is also interesting to compare our results with results for d = 4. For d = 4, axisymmetry is much 
less restrictive than for d > 4. This is because, the action of a 1-dimensional group such as SO (2) 
need not be orthogonally transitive, so the orthogonal decomposition ^ is not always possible, 
e.g. it does not apply for the Kerr solution. The natural d = 4 analogues of d > 4 axisymmetric 
spacetimes are spacetimes with a spacelike, hypersurface-orthogonal Killing vector field. Coordinates 
can then be chosen so that the Killing field is d/d(j) and there is a discrete isometry (p —cj). The 
metric then can be written in the form with dil^ = dcjP'. So this is the class oi d = 4 spacetimes 
analagous to our spacetimes. Algebraically special d = 4 vacuum solutions (with A = 0) with these 
symmetries were first classified by Kramer and Neugebauer [30], so we shall refer to them as KN 
solutions. 

By the Goldberg-Sachs theorem, the null congruence tangent to the repeated PND is geodesic 
and shear-free. For KN solutions it can be shown that it is also rotation- free, i.e., hypersurface- 
orthogonal [30]. Hence a KN solution belongs to the Robinson- Trautman (RT) or Kundt family of 
solutions depending on whether the congruence associated to the repeated PND is expanding or not. 
In 4d, the general vacuum solution belonging to either of these classes involves arbitrary functions of 

^ This answers a question posed in [1^: do there exist d = 5, A = 0, type D solutions with non-geodesic multiple 
WAND? For the type D examples with non-geodesic multiple WANDs encountered in our analysis, the WANDs are 
not discrete so they are not PNDs. So one could restate the question as; does there exist a d = 5 vacuum type D 
spacetime with a non-geodesic repeated PND? 

®We are taking results from the arXiv version of Ref. [5 , which differs significantly from the published version. De 
Smet worked in Euclidean signature and hence could not distinguish between the Schwarzschild solution and a KK 
bubble. 
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CoiOj = 











type II 


Coioj = 


Coijk 


= 






type D 


CoiOj = 


Coijk 


= Ciiij 


= Clijk = 





type III 


CoiOj = 


Coijk 


= Cijkl 


= Coiij = 





type N 


CoiOj = 


Coijk 


= Cijkl 


= Coiij = 


Clijk 


type 


Cal3"fS = 


= 









Table 1: Conditions for the various algebraic types. Note that the definition of type D involves 
secondary classification, i.e., consideration of n as well as i. 

time, and cannot be written down in closed form [6]. However, with the KN symmetries, the general 
Kundt (but not RT) solution can be obtained in closed form [30]. Special cases of RT solutions with 
the KN symmetries are the A-metrics and the C-metric. The Kundt family of solutions contains the 
B-metricsllI 

The main differences betwen these d = A results and our results are (i) the absence of time- 
dependent axisymmetric d > 4 RT solutions; (ii) the absence of a d > 4 analogue of the C-metric. 
The first difference extends beyond axisymmetry: the general class of d > 4 RT solutions was 
investigated in Ref. |T5] and found to be considerably simpler than the d = 4 class. In particular, 
the only d > 4 RT solutions with non-vanishing "mass function" are simple static generalizations 
of the Schwarzschild solution, in contrast with the d = 4 case where such solutions are generically 
time-dependent. Concerning point (ii), to explain what we mean by "analogue", we note that the 
main interest in constructing such a solution is to obtain an exact solution describing a black hole on 
a Randall-Sundrum brane, as explained in Ref. [31]. Such a solution would have d = 5, A < and 
would be axisymmetric (if the black hole were spherically symmetric on the brane) and describe an 
object with an event horizon accelerating along the axis of symmetry. The d = 4 C-metric belongs 
to the Weyl class, the RT class and the class of type D metrics. However, no d > 4 analogue was 
found in the generalized Weyl class (for A = 0) [4| or, as we have just discussed, the d > 4 RT class 
[Is] . Our results demonstrate that this negative conclusion extends to the d > 4 type D class too. 
However, we note that a type D metric of the form ([5]) does describe the near-horizon geometry of 
an extremal charged Randall-Sundrum black hole [24\ . 

1.3 Algebraic classification in higher dimensions 



For convenience, we shall review briefly the CMPP classification scheme [9]. This involves a null 
basis Co = i, ei = n, = mj, i = 2 . . . d — 1 where £ and n are null vectors obeying i - n = 1, and 
are an orthonormal set of spacelike vectors orthogonal to i and n. Consider a new frame related by 
a boost 

£ = X£, h = \~^n. rhi = mi. (6) 

A covariant tensor component Ta^^^ ap is said to be of boost weight s if its value in the new basis is 
related to its value in the old basis by 

Tai...dp = A Ta^ a^. (7) 

We are primarily interested in the Weyl tensor. The components of highest boost weight are the 
components Coioj, which have s = 2. The null direction £ is called a Weyl aligned null direction 

^We have ignored a special case that arises in the KN analysis, which occurs when the repeated PND is not invariant 
under the discrete symmetry (and must therefore be mapped to another repeated PND so the spacetime is type D). 
KN showed that the only such solution is the k = B-metric. However, this admits a second hypersurface-orthogonal 
spacelike Killing field, and the associated discrete symmetry does preserve the repeated PNDs. This implies that the 
solution also belongs to the Kundt class. 
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(WAND) if, and only if, these components vanish. This is independent of how the other vectors of 
the basis are chosen. For d = 4, WANDs are the same as PNDs. However, for d > 4, in general, no 
WAND exists and the spacetime is called type G. If a WAND does exist then the solution is called 
algebraically special. Another important difference in higher dimensions is that WANDs need not 
be discrete. The terminology PND is reserved for the case in which there are finitely many WANDs. 

If I satisfies the condition that all Weyl tensor components of boost weight 2 and 1 vanish then I 
is a multiple WAND. Using the tracefree property of the Weyl tensor, the multiple WAND condition 
is Coioj = Coijk = 0. If an algebraically special spacetime does not admit a multiple WAND then 
it is called type I. If it does admit a multiple WAND then it is type II, or perhaps more special. 
The more special types are: type III if all Weyl components of boost weight 2, 1, vanish, type N 
if all components of boost weight 2, 1,0, —1 vanish and type O if the Weyl tensor vanishes. The 
explicit conditions for these more special algebraic types are given in table [H The algebraic types 
are mutually exclusive, i.e., a spacetime of type III is not also type II. 

So far, we have discussed only the WAND which gives the so-called primary classification. One 
can then perform a more refined secondary classification by examining whether, or not, for £ given 
by the primary classification, it is possible to choose n to make further Weyl tensor components 
vanish. For example, if spacetime is type I with WAND d. and it is possible to choose n so that 
Ciiij = then the spacetime is said to be type Ij. We shall not make use of secondary classification 
except in defining type D spacetimes. A spacetime of primary type II, with multiple WAND is 
said to be type D if n can be chosen so that Cmj = Cujk = 0. Note that this implies that n also is 
a multiple WAND. 

We are interested in solutions of the vacuum Einstein equation ([1]). For such spacetimes, the 
WAND conditions can be reformulated in terms of the Riemann tensor: 



Rq^qj = 0^ WAND, Roioj = Roijk = <^ multiple WAND. (8) 

We shall make use of several general results for warped product spacetimes. A warped product is a 
spacetime of the form 

ds^ = A{yfgab{x)dx^dx' + B{xfg,,{y)dy'dy^, (9) 

where gab is Lorentzian and gij is Riemannian. Such a spacetime is type D or O if the Lorentzian 
factor is (i) two-dimensional; (ii) a three-dimensional Einstein space; (iii) a type D Einstein space 



2 Static, axisymmetric, solutions 

In this section we consider higher-dimensional solutions that are static and axisymmetric, i.e., they 
admit a hypersurface orthogonal timelike Killing vector field that commutes with the generators of 
SO{d — 2). Introduce coordinates adapted to the isometries: 

ds^ = -A{r, zfdt^ + B(r, zfidr"^ + dz^) + C(r, zfdn'^. (10) 

The components of the Einstein equation ([T|) for this metric are given in Ref. [32]. Define a complex 
coordinate w = {r + z2:)/\/2. Consider R^. — R^ + 2iRZ, = 0. This gives 

d'^A dAdB fd'^C dBdC\ , , 

where d = d/dw. This implies 

dB _ C'^-^d'^A + {d- ?,)AC'^-^d'^C 
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We must consider the denominator since it could vanish identically, i.e., AC^ ^ might be constant. 
The equation Rj + {d - 3)Rgl = {d - 2)A implieiil 

^g,_3 = [(^ - - 4) - (d - 2)AC2] . (13) 

Constancy of d{AC'^~^) implies the RHS must vanish hence C is a constant and A is positive. But 
then A must also be constant so the spacetime has a fiat time direction. This is incompatible with 
positive A. Hence AC^'^ cannot be constant. 

If A and C are known then equation (jl2p determines B. Furthermore, i?^ + i?f — 2R\ = implies 

{-^Bf _AB_ {d-3)AC _AA_ ^. VA-VC 
52 B ~ 2C 2A AC ■ ^ ' 

It can be checked that this is compatible with equation (fT2|) . 

We assume that the spacetime is algebraically special, so it admits a WAND i. We shall assume 
for now that the WAND is axisymmetric. Assuming d > 4, this implies that it is orthogonal to S'^~^. 
By rescaling i we can arrange that 

n d Af , ^9 . , 

dt^B ""'d^^ s^'^f^yr, z)—y (15) 

for some function a{r,z). Staticity implies that 

d A f d d \ 

""^di'B (^^"^"('■'^^a;: +^™"('''^)a^J ' (1^) 

is also a WANe[^, i.e., WANDs come in pairs, which implies that the algebraic type must be I, D or 
O [H]. Choose 

= ^-sinQ(r,2;)^ + cos a(r, z)^^ , (17) 

and 

m°' = ^e", a = 3,...,d-l (18) 
where e" is a vielbein for S"^^^. We find that the WAND condition ([8|) reduces to 

Re (e^^^VF) = X, (19) 

where 

d^A d'^C dAdB BBdC ,^ C^-^{d''AdC-dAd''C) 

^ _ Ai? , AC • VC , {VBf _ {d- 2)A^ ^ VC 

2A B ^ 2C AC + 52 - 2C ^ ' ^ ^2 A 
where the second equality in each case follows from equations ()12p and ()14|) . The spacetime is 
algebraically special if, and only if, there exists a real solution a of the WAND condition. Hence 

> |X| 44> algebraically special (22) 

Now consider the additional condition required for a multiple WAND (equation ([5])). This gives the 
single equation Im (e^*"W^) = 0. Combining with the type I condition gives 

e^^'^W = X. (23) 



* We have defined V = (dr,dz), A — V'^, and indices are raised with the flat metric dr^ + dz^ . 

^ These nuU vectors don't obey £ ■ n = 1 but this can arranged by rescaling them, which doesn't affect anything 



below. 
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We conclude that, assuming an axisymmetric WAND 



\W\ < \X\ ^ TypeG 

\W\ > \X\ ^ Type I (24) 
\W\ = \X\ ^ Type D or O 

We shall now comment on our assumption that the WAND is axisymmetric. In general, this need 
not be true. However, for odd d, is even dimensional so the projection of £ onto the sphere 

must vanish somewhere. Working at such a point, we can argue as above to arrive at equations ()24|) 
that depend only on r and z and must therefore hold everywhere on the sphere, which implies the 
existence of an axisymmetric WAND. Hence, for odd d, there is no loss of generality in restricting to 
an axisymmetric WAND. For even d, this argument does not work. However, in section [U we shall 
consider axisymmetric spacetimes with a non-axisymmetric multiple WAND, and show that no such 
spacetime is static and axisymmetric0 Therefore the multiple WAND of a static axisymmetric type 
D spacetime must be axisymmetric. However, it is possible that some spacetimes with even d > 4 
and \W\ < \X\ may be type I with a non-axisymmetric WAND. 

To illustrate these conditions, consider the z-independent d = 5, A = solution of Ref. [33] . 
written in the form given in Ref. |34j 

l + m/Rj \l + m/R) 



where a = y^/?^ + /3 + 1. Assume m 7^ (so the spacetime is not flat). Then a calculation reveals 
that \W\ = \X\ if, and only if, /? = or /? = 1. The first possibility gives the Schwarzschild black 
string. The second possibility gives a boost invariant singular spacetime discussed in Ref. [21j . This 
spacetime is of the form ([5]) discussed in the introduction. 

Another interesting example is the static Kaluza-Klein bubble (the product of a flat time direction 
with the Euclidean Schwarzschild solution). This can be obtained by taking the limit /? — > cx) of the 
above metric. This spacetime has = 0, X 7^ hence it is type G. 

Since type G is distinguished from type I only by an inequality, it is possible that there exist 
(connected) analytic spacetimes that are type G in some open subset of spacetime and type I in 
some other open subset. Indeed, if we choose m > and /? = 1/2 in the above metric then it is type 
G for i? ~ m but type I for ^ m. As discussed in the introduction, this kind of behaviour 
suggests that the type I condition alone will not be much help in solving the Einstein equation. 

It would be nice to use the type D condition obtained above to solve the Einstein equation. 
However, we have not made progress using the coordinates employed here. (Even in d = 4, this 
approach would not work for A 7^ 0.) However, in subsequent sections we shall see that all static 
axisymmetric type D solutions can be found, indeed we shall relax the condition of stationarity and 
determine all axisymmetric type D solutions. 



3 Axisymmetric solutions with an axisymmetric WAND 
3.1 Introducing coordinates 

In this section we shall consider general (possibly time-dependent) axisymmetric spacetimes with 
an axisymmetric WAND. First we shall introduce coordinates adapted to the WAND. Consider the 

'^'^ Actually, the spacetimes we find there are "static" and "axisymmetric" but they are not "static and axisymmetric" 
because the generator of time translations does not commute with the generators of axisymmetry. 

^^The general behaviour appears to be that, for < /3 < 1 and m > 0, the solution is type G near R = m and type 
I for m. For (finite) /3 > 1 and m > 0, it is type I near R = m and type G for R^ m. 
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metric in the form where spacetime is locally a warped product M3 x S'^ ^. Axisymmetry implies 
that the WAND is tangential to M3. 

We shall choose the local coordinates on M3 as follows. Pick a 2-surface in M3 transverse to 
the WAND I and let be coordinates on this surface, where a = 1, 2. Now carry these coordinates 
to the rest of spacetime along the integral curves of I, and let f be the parameter distance along 
these curves. Now use {x"',f) as coordinates on M3, so ^ = 8/ dr. The metric takes the form 

ds'^ = 2gradrdx^ + g-^^dx^dx^ + E'^dn'^, (26) 

where g^a 7^ for some d. Without loss of generality we may assume gfi 7^ 0. Now let r = 
/ df li'f^, x"')df, V = x^, z = x^ and use coordinates {v,r,z). In this chart, the metric takes the form 



dr — -U {v, r, z) {dv + B{v, r, z)dz) + C{v, r, z)dz 



ds^ = 2[dv + B{v,r,z)dz] 

+ D{v,r,zfdz'^ + E{v,r,zfd^\ (27) 



for some functions U,B,C,D,E. The WAND is proportional to d/dr so we can rescale it so that 
i = d/dr. It is convenient to complete this to a null basis as follows: 

Ldx'^ = dv + Bdz, n^dx^ = dr-'-ULdx^ + Cdz, 

= Ddz, = Ee", (28) 

where e" (a = 3 ... d — 1) is an orthonormal basis of 1-forms on S'^~^ with no dependence on v,r,z. 
We shall denote the spacelike basis 1-forms collectively as m*, i = 2 . . . d — 1. 

Now consider the null congruence associated with the WAND. This is geodesic if, and only if, 
drB = 0. The "expansion matrix" of the congruence is 

Sij = mfm^- V(/^) = diag ( 1 . (29) 

The expansion of the congruence is the trace of this matrix and the shear tensor is the traceless 
part. The rotation matrix of the congruence vanishes: this is a consequence of axisymmetry. 

3.2 Type II or D implies geodesic or dS^^ x 5''^"^ 

Consider the case that £ is a multiple WAND. It has been shown that for a type III or N Einstein 
spacetime, the multiple WAND is geodesicQ [11]. In this section we shall prove a generalization of 
this result for axisymmetric vacuum spacetimes with an axisymmetric multiple WAND, of principal 
type II, i.e., the spacetime is type II or D. 

Assume that spacetime is type II or D and that the multiple WAND is axisymmetric and not 
geodesic. Define the matrix 

^^J=CQ^lj. (30) 

Let # = Define also 

Li = m'ttVJ^, (31) 

so the WAND is geodesic if, and only if, Li = 0. It has been shown [12J that if the multiple 
WAND is not geodesic then ^>jj must be symmetric with an eigenvalue equal to — with associated 
eigenvector Lj. In our case, axisymmetry implies that L^ = 0. Hence we must have ^22 = Now 
axisymmetry implies that = and oc 5^^, so we must have 

^ij = $ diag (-1, 2/{d - 3), ... , 2/{d - 3)) . (32) 



^See footnoted 
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For any type D spacetime satisfying ([T]), equation (27) of Ref. [12\ relates to Sij: 

= -5'fcfc$(ij) + ^Sij + ^(ik)Skj + ^(jk)Ski + 3 {^[ik]Skj + ^[jk]Ski) + CikjiSik (33) 

It has recently been observed that this equation is actually valid for any solution admitting a multiple 
WAND, i.e., it is also valid for type II [35j. Using the above expressions for <!>,-,• and Sij, this equation 
reduces to 

= 0. (34) 

If $ = then the traceless property of the Weyl tensor implies that Cikjk = 0. This implies that 
C2a2a = Ca^y/Bj + C2a2f3 = 0. But, by axisymmetry, C2a2/3 oc dap and Cap-fS oc {5a^5ps - Sp^Sas), so 
these equations imply that C2a2/3 = Cafj-yS = 0, i.e., Cijki = 0. Therefore all Weyl components of 
boost weight zero vanish, implying that the spacetime is type III (or more special), a contradiction. 
Hence <I> 7^ 0, so we must have 

drE = 0. (35) 

We now substitute this result into the equation RoaOis = following from the WAND condition ([8]). 
This gives drB {d^E - Bd^E) = 0. But drB / (the WAND is non-geodesic) hence 

d^E - Bd^E = 0. (36) 

Taking a r-derivative of this we obtain d^E = 0. Plugging this back into the equation gives dzE = 0. 
Hence E is constant. Therefore the spacetime is a direct product M3 x S'^~^. This is only possible if 
A > 0, then the Einstein equation implies that the solution must be dS^ x S'^~^ (which is type D). 

In summary, we have shown that an axisymmetric multiple WAND of an axisymmetric type II 
or D vacuum spacetime must be geodesic unless the spacetime is dS^ x S'^^^ . Combining this with 
the results for type III or N [11], we learn that an axisymmetric multiple WAND of an axisymmetric 
vacuum spacetime must be geodesic unless the spacetime is dSs x S'^~^ or type O. 



3.3 Solutions with a geodesic WAND 

The results Ref. [llj and the previous subsection establish that an axisymmetric multiple WAND 
in a vacuum solution is always geodesic (unless the spacetime is dSs x S'^~^ or type O). In this 
subsection we shall determine all solutions with an axisymmetric geodesic WAND. We shall not 
assume that the WAND is a multiple WAND (so a priori the solution might be type I but we shall 
see that this does not happen). 

We have drB = because the WAND is geodesic. We can now introduce new coordinates v' and 
r' such that v = v' + F{v',z), r = r'G{v',z) for some functions F,G that can be chosen to bring 
the metric to the same form as before but with B = 0. Dropping the primes on the coordinates, we 
have 

ds^ = -U{v, r, z)dv^ + 2dvdr + 2C{v, r, z)dvdz + D{v, r, zfdz^ + E{v, r, zfdfl'^. (37) 

By rescaling, the WAND can be taken to be ^ = d/dr. We saw above that the null congruence 
associated with the WAND has vanishing rotation. Since it is geodesic, this implies that it is 
hypersurface orthogonal. In the above coordinates, it is orthogonal to hypersurfaces of constant v. 
Furthermore, r is an affine parameter along the null geodesies!^ There is some coordinate freedom 
remaining: the form of the metric is invariant under the transformations 

V — > y{v)., r r/dyV, 

r r-F{v,z), (38) 

z — > z{v,z). 

^^If we assumed that ^ is a multiple WAND then the r-dependence of the metric could be read off from Ref. [14] , 
However we shall not make this assumption. 
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All of this is well-known in the context of 4d solutions with a hypersurface orthogonal null geodesic 
congruence |6j. 

We shall employ the same null basis as before (i.e., ()28p with B = 0). The Riemann tensor of 
the above metric in this basis is given in Appendix |Al The WAND condition ^ reduces to 

d'^D = dlE = 0, (39) 

Hence 

for some functions Dq, Di,Eq, Ei. The 00 component of the Einstein equation is now automatically 
satisfied. Axisymmetry implies that the Oa component is trivial. The 02 component reduces to an 
equation linear in C: 

drC-[^^-id-3)^jdrC- 2{d - 3)-^C = 2{d - 3) ^— ^ J . (41) 

This will determine the r-dependence of C. There are several different cases to consider. 
3.3.1 ^1 / 0, L>i = 

We can use the residual freedom in r and z (equation (f38l) ) to set £^0 = and = 1, i.e. D = 1. 
Then ([H]) reduces to 

5.^C + ^a.C = 2^^Ml, (42) 
r r El 

which can be solved to give 

Civ, r, z) = Coiv, z) + + -^r, (43) 

for arbitrary functions Cq and Ci. The r-dependence of U{v, r, z) is determined by the 01 component 
of Einstein's equation: 

where Uq and C/i are arbitrary functions, and 

{log(r) , d = 5 

1 • (45) 

(d-5)r'^-5' ""^^ 

The r-dependence of the metric is now fully determined. Comparing coefficients of terms with 
different r dependence in the remaining components of the Einstein equation can be used to restrict 
the arbitrary functions above. The a/3 components of Einstein's equation give 

Ci = 0, a,Co = 0, (46) 

C^o(^,^) = ;^^(^-C^o(^)'> (47) 
The residual coordinate freedom z ^ z — f{v) can be used to set 

Co = 0. (48) 
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The 22 component of Einstein's equation gives 

El El d ■ 

Now the 12 component of Einstein's equation imphes 



A 



0. 



a^C/i + 2((i-3)^^C/i = 0, 



which in turn imphes 



Ui{v,z) 



El 
m{v) 



Ei{v,zf('i-^y 

for some arbitrary function m{v). The 11 component of Einstein's equation then imphes 

d^Ei{v,z) 



dy'm{v) = {d — A)m{v)- 



EUv,z) 



(49) 

(50) 
(51) 

(52) 



m{v) = imphes that the spacetime is conformahy flat (type O), so assume that m{v) ^ 0. Then 
([^ imphes 

\m{v)\y(^-^) 

Ei{v,z) = , 53) 

for some positive function g{z). Inserting this into equation (j49p gives the hnear equation 



g"{z) + j^g{z) = 0, 



where L is defined by 



A 



id-l)e 
L2 ' 



eG {-1,0,1}. 



Define a positive constant n via the first integral 



g'izf + ^g{zf 



where ?7G{— 1,0,1}. Using the freedom to shift z by a constant (and z 
where G{z) is given by 



(54) 
(55) 
(56) 

-z) we have g{z) = iJ.G{z), 



G{z) 


1] = \ 


r? = 


rj = —1 


e = 1 


L sin(z/L) 






e = 


z 


a 




e = -1 


L sinh(z/L) 




L cosh.{z/L) 



where a is a positive constant. Defining new coordinates (V, p) by 



dV = fi- 



dv 



the metric becomes 



ds' 



dz"^ + G{z)'' 



m(?;)|V(rf-4)' 
M 



1 



4-4 



r|m(i;)|i/('^-4) 



rjp^ dV^ + 2dVdp + p^dn^ 



(57) 
(58) 



where M = sign{m{v)) p'^~^ . This is the warped product of a line, parametrized by z, with the 
(c? — l)-dimensional Schwarzschild (anti)-de Sitter metric. The e = rj = case is the Schwarzschild- 
Tangherlini black string. The e = 0, r] = 1 case corresponds to taking a de Sitter slicing of Minkowski 
spacetime and replacing the de Sitter slices with the Schwarzschild-de Sitter metric. The e = rj = 1 
case corresponds to doing the same thing for de Sitter spacetime. The e = — 1 cases correspond to 
the same idea for slicings of anti-de Sitter space (In d = 5, e = —1, rj = is the AdS black string 
of Ref. [E]). A warped product whose Lorentzian factor is a type D Einstein space is also type D 
[12) . Hence these solutions are all type D. 
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3.3.2 El / 0, Di^ 0, Dq = 0: Robinson-Trautman solutions 



The coordinate freedom (l38|l can be used to set £"0 = 0. From (j29[) . these solutions have vanish- 
ing shear and non-vanishing expansion. Therefore they belong to the class of higher-dimensional 
Robinson-Trautman solutions [15]. To give a self-contained presentation, we shall rederive these 
solutions here (with the additional restriction of axisymmetry) . Using a transformation z z'{z, v), 
we can set Di = 1. The general solution to equation (j4ip is 

C{v,r,z)=Ci{v,zy + ^^^, (59) 



r 

where Ci and C2 are arbitrary functions. But now R22 — Raa (no sum on a) is independent of U and, 
by the Einstein equation, must vanish. Equating coefficients of terms with different r-dependence 
gives C2 = and 

a.c, - Mic, + ^ = 0^ (61) 

The 22 component of the Einstein equation can then be solved to determine IJ: 

.) = - ^fl + 2.a,C, - {ci + ^) , (02) 

where C/i is an arbitrary function. Now examining the 12 component of the Einstein equation gives 
dzU\ = 0. The 11 component of the Einstein equation reduces to 

Eid^Ui + {d-l) {Uid^Ei - CiUid.Ei) = 0. (63) 

If C/i = then it can be shown that the above equations imply that the Weyl tensor vanishes hence 
the solution is type O. If C/i 7^ then we can use the gauge freedom v — > V{v) and r — > r/d^V to 
set Ui = ^ for some non-zero constant /x. Then ()63p gives d^Ei = CidzEi. From ()6ip we then learn 
that dzCi = 0. This implies that Ci can be gauged away by a shift z ^ z — f{v). In the new gauge 
we have d^Ei = 0. The solutions of (pOj) are Ei = Rsin{z/R), z or Rsmh{z/ R) (using z — > z — const 
and z —z to simplify) where i? is a positive constant. R can be set to one by rescaling z, v and 
r. The solution takes the final form 

ds^ = -(^k- - -^f^r^^ dv'^ + '^dvdr + r^dS^, (64) 

where M is a non-zero constant, and dS^ is the metric on a d — 2 dimensional space of unit constant 
curvature of sign k. This generalized Schwarzschild metric is of type D |15j . 



3.3.3 El 7^ 0, Di^{), Do^O 

We use the transformations ()38p to set Di = 1 and Eq = 0. The general solution to equation ()4ip is 
C{v,r,z) = {Doiv,z)+rf(^Ci{v,z) + C2{v,z) J ^.-3(^^ ^ ) " (65) 

where Ci and C2 are arbitrary functions. Now we consider the 22 and aa components of the Einstein 
equation. These equations are linear in U and drll and can be solved algebraically to determine U 
and drU. The r-dependence is completely determined hence consistency of the solutions for U and 
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drlJ gives an equation whose r dependence is fully determined. Equating coefficients of terms with 
the same r dependence then gives C2 = together with 



d.Ei = -J^DoEi - Eid.Ci + Cid.Ei. (66) 



The solution for U is then 



U{v,r,z) = -C^Do + 2CiDo— ^ 



E( ^ " ^ " ^1 Ef 
+ (2a,C7i-^Z)o-2C?i?o)r-(^ + Cf)r2. (67) 

These results imply that the Weyl tensor vanishes. Hence these solutions are type O, i.e., Minkowski 
or (anti)-de Sitter spacetime. 

3.3.4 El = 0: Kundt solutions 

Solving ([^T]) gives 

dzEo\ , 1^ , , , 2 



C{v, r, z) = Co{v, z) + [C,{v, z)Do{v, z) + 2{d - 3)^ j r + -C,{v, z)Di{v, z)r\ (68) 

where Cq{v,z) and C^{v,z) are arbitrary functions. The aa component of the Einstein equation 
does not involve U so its r dependence in completely determined. Equating coefficients of terms with 
different dependence on r gives D\[{d — ^) — KEq] = 0. Hence either Di = QotEq = (d — 4)/A. The 
latter implies that spacetime is a direct product M3 x S"^~^, which requires A > and the spacetime 
must then be locally (iS'3 x S'^~^ which is of algebraic type D. 

Assume instead that Di = 0. We now have drD = drE = so the geodesic congruence is free 
of expansion and shear as well as twist. Spacetimes with vanishing expansion, shear and twist are 
referred to as Kundt spacetimes [6l|36]. All vacuum Kundt solutions are type 11 or more special for 
any d > 4 [T7]. General d-dimensional Kundt spacetimes have been discussed recently [37]. The 
general solution cannot be obtained in closed form. We shall now analyze such solutions assuming 
axisymmetry, which enables further progress to be made. 

With Di = 0, we can use the transformation z z'{v, z) to set Dq = 1. Write the solution for C 
as C{v, r, z) = Cq(v, z) + rCx(v, z). The shift r^r- F{v, z) has the effect Cq ^ - d^F - CiF. 
Hence we can choose F{v, z) to set Co = 0. To summarize, we have brought the metric to the form 

ds^ = -U{v, r, z)dv^ + 2dvdr + 2rCi{v, z)dvdz + dz^ + Ea{v, zfdfl^. (69) 



Some gauge freedom remains. The transfomations of the form (I38p that preserve this form of the 
metric are 

v = v', z = z' + f{v'), r = r' + g{v',z'), d,,g + Cig = -8,, f , (70) 

v = V{v'), r = jf—, z = z'. (71) 
The aa component of the Einstein equation reduces to 

^ + «<-4»-C.M2-y+A^0. (72) 
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and the 22 component of the Einstein equation reduces to 



- icf - (d - 3)^ - A = 0. (73) 

The 01 component of the Einstein equation is satisfied if, and only if, 

U(v,r,z) = \J^(v,z)^TlJ^(v,z)^T''U2{v,z), (74) 
where f/o ai^d ^\ arbitrary functions, and 

^2{v.z) = \d.C, + ^Ci^ - \cl - A. (75) 

The 12 Einstein equation reduces to 

d,U, = -d,C, -{d- 3)Ci^ - 2id - 3)^4^. (76) 

Finally, using the above equations, the 11 Einstein equation reduces to 

dPo +(c, + {d-3)^) d.Uo + ( d.C, + {d- 3)C,^) C/o = (d - 3) f 2^ - U, ^"^^ 



Eq J \ Eq J \ Eq Eq 

(77) 

As is famihar for Kundt solutions, the equations of motion separate into the "background" equations 
()72p and (j73p . which must be solved to determine Eq and Ci. Given a solution of these equations, 
the other equations can be integrated to determine Uq and Ui. The second step is trivial because the 
equations are linear. Hence solving the background equations is the non-trivial step that remains. 
However, the general solution to the background equations is not known analytically. 

Since the background equations do not involve f -derivatives, solving them equations is equivalent 
to solving the corresponding equations assuming that Eq and Ci are independent of v and Uq = 
Ui = 0. But in this case, the metric is static. In fact, we shall see below that the general type D 
axisymmetric Kundt metric is of this form. The background equations can only be solved in special 
cases e.g. the general solution with d = 4 can be determined, and the general solution with d = 5, 
A = and U2 = can also be obtained [21j. Some time-dependent solutions based on the latter 
solution of the background equations were obtained in Refs [26l [27] . 

It is convenient to define a positive function W{v, z) by 

W{v,z)=Wo{v)eM- f Ci{v,z')dz'), (78) 
where Wq{v) is an arbitrary positive function, so 

The background equations become 

dlE, {d^Eof d.Wd.Eo d-A 

W'^"^-^^^-^-^- 



Equation ([77|) becomes 



dlW {d-2,)d,Wd,EQ 

^^ = -w WW, ^- 
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These equations imply that 
hence we can always choose Wo{v) so that 

U2{V,Z) 



dz ( W2) = 0, 



with k £ {1,0,— 1}. We can now define a new coordinate R by 

r = W{v,z)R. 

The metric becomes 



(83) 



(84) 



(85) 



ds^ = W{v,z) 



Uo{v,z) 
W{v,z) 



+ R[ Ui(v,z) 



W 



kR^ 



dv'^ + 2dvdR } +dz'^ + Eoiv, zfdn^ 



(86) 

We now consider a further classification of the Kundt solutions according to their algebraic type. 
Using the above equations to simplify the Weyl tensor, we find that the only independent nonzero 
components are: 

d,Eo A 



Cv 



2En 



d-l 



(87) 



C 



la/32 



Ca/375 — 26a[sS'y]f3 



Ci 



d- 1 



A ^{d^Eo^ 



K2 
^0 



E2 
-^0 



En 



+ 2- 



En 



2{d-3) 



c 



1 



2{d - 3) 



a,(a,[/o + CiC/o) + 



Eq \^ 



En 



(90) 



Note that while Cqioi, ^0112, ^02121 (^1212 and 6*^2/32 are nonzero, they are related to the above 
components by the tracefree property of the Weyl tensor. The first two Weyl components written 
above are of boost weight 0, while the remaining two are of boost weight -1 and -2 respectively. 
Hence the solutions we are considering here are at least type II, confirming the general result of Ref. 

Type III and N 

Consider the case in which the solution is type III, or more special. In this case, the Weyl 
components of boost weight vanish. This gives, for d > 4, the following equations 



C 



dzEp 

Ea 



{d,Eo) 
El 



2A 



A 



1 



El {d - 1) 



(91) 
(92) 



Note that equation (|92p is not present for d = 4. Solving this equation gives 

'Lsin(z/L) ifA>0, 
Eo{v,z) = Iz if A = 0, 

^Lsinh(z/L) if A < 0, 



(93) 
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where L > is defined by (I55p , and we have used the freedom (I70p to ehminate an arbitrary function 
of V (we've also fixed signs using z zizz). Equation ()9ip now determines Ci: 



Ci{v,z) 



(2/L)tan(z/L) if A > 0, 

if A = 0, (94) 

^-{2/L)tanh{z/L) if A < 0. 



Since Eq is independent of v, equation (f89l) gives that Cia/32 = 0, hence these solutions are type N 
or O. There are no axisymmetric type III Kundt solutions for d > 4. However, such solutions do 
exist for d = 4 [6] . 

Continuing the analysis, note that the coefficient of r in Ciaip (given by ([90]) ) vanishes, and so 
Ciai/3 reduces to 

Cic.w = ~ 2(d-3) ^^^^"^° + ^^^^ 



U2{v,z) can be calculated using ([75|) : 

U2{V,Z) = 



'-(1/L2)sec2(z/L) ifA>0, 

if A = 0, (96) 

_(1/L2)sech2(z/L) ifA<0. 



Equation ([76]) implies that J7i(t;,z) is independent of z: Ui{v,z) = Ui{v). We can then use the 
transformation v y{v), r — > r/d^V to arrange that 

Ui = 0. (97) 

Finally, equation ()77p can be solved to determine Uq. For A = 0, the solution is, for d> A 

F(v) 

Uo = ^ + G{v), (98) 

(for d = 4, the first term is replaced by F{v) log z). A shift r ^ r — /(i;) can be used to set G{v) = 0. 
The only independent non-zero component of the Weyl tensor is (|95p . This reveals that the solution 
is type O if, and only if, F{v) = 0. Therefore, the general axisymmetric type N Kundt solution with 
A = 0, is given by the following metric (for d > 4) 

ds^ = -^^dv"^ + 2dvdr + dz'^ + z^dn"^. (99) 

z"- * 

The null vector field i is covariantly constant, and so the solution above belongs to the family of 
pp- waves [38] . 

For A < 0, the solution for Uq is 

Uo = cosh2(z/L) [F{v)I-{z) + G{v)] , (100) 

where 

^- -2 = / 7T^^ ' 101 

Define a new coordinate R by 

r = Rcosli^iz/L). (102) 

The metric becomes 

ds'^ = cosh^iz/L) [-{F{v)I^{z) + G{v) + R^/L^)dv'^ + 2dvdR] + dz'^ + smh^ [z / L)d^^ . (103) 
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Now the transformations R R — f{v) followed hy v ^ ^ ~^ R/d^V can be used to eliminate 

G{v), giving the final form of the solution: 

ds^ = cosh^{z/L) [-{F{v)I-{z) + R^ / L^)dv'^ + 2dvdR] + dz^ + smh^{z/ L)d^^ . (104) 

A similar analysis for A > (or L iL) gives 

ds^ = cos^{z/L) [-{F{v)I+{z) - R^/L^)dv^ + 2dvdR] + dz"^ + sm^{z/ L)dn^ , (105) 

where 



dz 

cos2(z/L)"lh^'^-3(2/L)' 

These solutions are type O, i.e., isometric to (anti-)de Sitter space, if, and only if, F{v) = 0. If F{v) 
is not identically zero then these metrics are the general axisymmetric type N Kundt solutions for 
A / 0, d > 40 It seems natural to interpret the type N solutions as describing gravitational waves 
propagating in a type O background. 

Type D 

Now consider type D solutions, for which there exists a second multiple WAND n' . Note that 
n' need not coincide with the null basis vector n defined above. If n' were not axisymmetric then 
the solution would be encompassed by the analysis of section [H However, the results of that section 
reveal that, in this case, both multiple WANDs would fail to be axisymmetric, which is not the case 
here. Hence we can assume that n' is axisymmetric. The most general form it can take is 

n' = n — 2^(^' ^1 -^)^ ^ + '^(^' ^' ^) "^2 (10''') 

where a{v,r,z) is arbitrary. Let us change to a new null frame {i',n',m^), with 

£' = i, m'2 = 1712 — a£, m'^ = rua- (108) 

Note that a = corresponds to the frame used above. The fact that £ is a multiple WAND guarantees 
that Weyl components of boost weight are the same in the two frames. The negative boost weight 
components in the new frame are related to the components in the old frame by 

^laip ~ C'l«l/3 — 2aCiQ,^2 — + 0?C2a2p- (109) 

Type D solutions are those for which Weyl tensor components of boost weight —2 and —1 (in the 
new frame) vanish, giving 

Cid^Eo + 2d,d,Eo + a {Cid.Eo + 2dlEo) = 0, (110) 

Clal/3 — 2aCiQ,/32 — + a^C2a2l3 = 0, (HI) 

where the second equation has not been written explicitly for brevity. Note that equation (jllOp 
implies either a = a(f,z); or C2a2(5 = C'ooi/? and = 0. In the latter case, equation (jllip 

implies Ci^ip = 0, and then one finds that n is a multiple WAND, i.e., one can set a = 0. Hence, 
in either case, we have a = a{v,z). Therefore, in equation (jllip . the only term with r dependence 
is that contained in Ciaip, which must vanish, giving (from equation (I90p ) 

2U2^ + {d,U, + 2d^C,) ^ = 0. (112) 



'For d = 4, these solutions are a special case of more general type N Kundt solutions discussed in Ref. [39j . 
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To simplify the analysis, assume that the spacetime is not dS^ x S'^~^ (which we already know is 
type D). The results of section 13.21 imply that the second multiple WAND n' must be geodesic. 
Axisymmetry implies that the geodesic equation reduces to 

m'^^n'^bn'a = 0. (113) 

The LHS is linear in r so this gives two equations: 

2ad^a + 2d^a + c?Cx - aUi + d^Uo + CiUq = 0, (114) 

a^^C/i + 2dyCi - 2aU2 = 0. (115) 
Proceed by simplifying equation (jll2p using equation (1115^ : 

«(„. 11)^0. (UO) 

Note that dzE^ is not identically zero, since otherwise eq. (j72|) implies that Eq is constant, which 
gives dS'i X S"^-3. There are two cases. 

Case 1. dyEo = —adzE^. Using this to eliminate d^E^ from eq. (jllOp gives dza = 0, so a = a{v). 
Now we can use a transformation of the form (j70p to reach a gauge in which dyE^ = 0, i.e., a = 0. 
Substituting ([73]) into (jllSp then gives dyCi = 0. Now equation ([76P gives dzUi = 0. Define a 
positive function W{z) by equation ([79]) . Equation ([75]) reveals that U2 is independent of v so 
equation ([83]) implies ([84]) as before (using the freedom to rescale by a constant). Equation ([114p 
gives Uq = F{v)W{z). Defining the coordinate R by ([85]) . the metric can be brought to the form 
([M]): 

ds^ = W{z) {- [F{v) + RUi{v) - kR^] dv^ + 2dvdR] + dz'^ + Eo{zfdn^. (117) 

The transformation R ^ R — G{v) can be used to set F{v) = 0, then a transformation v y{v), 
R — > R/dyV can be used to set Ui = 0. The metric is then 

ds"^ = W{z)dJ:'^ + dz^ + Eo{zfdn^, (118) 

where dH"^ is the metric on a 2d Lorentzian space with Ricci scalar 2k, i.e., Minkowski or (anti-) de 
Sitter spacetime. 

Case 2. U2 = 0. Equations ([75]) and ([73]) imply 



Eq Eq 



while equations (|115p and ([76]) imply 



a„(|ijj=2(d-3)^ (120) 



The integrability condition for these equations is 

hence d^E^ = h{v)dzEo for some function h{v). A gauge transformation of the form ([70p can be 
used to reach a gauge in which B^Eq = 0, i.e., h{v) = 0. Equation ([120p now gives 9t,Ci = and 
([76]) gives dzUi = 0. Now if a = then we are back to case 1, so assume a / 0. Then the coefficient 
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of a in equation (jllOp must vanish. But this is the case discussed below equation (jllip . where n is 
a multiple WAND, so one can set o = after all, leading back to case 1. 

In summary, we have shown that, for a general type D axisymmetric Kundt metric, one can find a 
gauge in which Eq and Ci are independent of v and Uq = Ui = 0. The metric can be transformed to 
the form (|118p Conversely, the warped product structure of (IllSI) implies that any such solution 
is type D or O Type O corresponds to the solutions for Eq{z) and Ci{z) found in our discussion 
of type N, i.e., equations ([93]), ([Mil . 



4 Axisymmetric solutions with a non-axisymmetric WAND 

Consider the Kaluza-Klein bubble spacetime [2^ (generalized to include A) obtained by analytic 
continuation of the Schwarzschild solution: 

ds^ = r^ds\dSd-2) + TTTT + U{r)dz\ U{t) = 1 - - (122) 
U[r) r" a — 1 

where m 7^ is a constant, and dSd^2 is {d — 2) dimensional de Sitter space: 

ds^{dSd-2) = -dt^ + cosh^ t dn^. (123) 

This spacetime is obviously axisymmetric. It is a warped product of dSd-2 and and is hence 
type D [12]. We did not discover this spacetime above. This is because the multiple WANDs live 
in the dSd-2 directions and hence must have non- vanishing components along S"^~^, i.e., they are 
not axisymmetric. Here this is possible because the axisymmetry SO{d — 2) is part of the bigger 
SO{l, d — 2) de Sitter symmetry. We shall show that this symmetry enhancement is necessary for a 
multiple WAND to be non-axisymmetric. 

Consider first the case in which we have a non-axisymmetric (multiple) WAND i that is every- 
where orthogonal to the S'^~^ orbits of SO{d — 2), i.e., the only non-zero components of the WAND 
(in the coordinates of are = i"'{x,Q), where refers to the coordinates on S'^~^. Now, since 
the Weyl tensor is axisymmetric, it is clear that i^{x,il,Q) is also a (multiple) WAND where 0,q is an 
arbitrary point on S'^~^. But this new WAND does not vary on S'^~^, i.e., it is axisymmetric. Hence 
we conclude that, if the WAND is everywhere orthogonal to S"^"^ then there is no loss of generality 
in assuming that it is axisymmetric. 

Assume instead that we have an axisymmetric spacetime with metric @ and that a WAND i is 
not orthogonal to S'^~^ at some point. Then the same must hold in a neighbourhood of that point. 
Consider the "unphysical" spacetime M3 x S"^~^ with the product metric obtained by multiping (|4|) 
by E{x)~'^. We shall work with this spacetime for most of this section. Obviously ^ is a WAND of 
this spacetime. By rescaling I we can ensure that the projection of I onto S'^"^ is a unit vector (in 
our neighbourhood). Hence we can write i = (eo + e^)/^/2 where cq is a timelike unit vector in M3 
and 63 a unit vector on S"^~^. Choose n = (— eo + e3)/\/2. Choose ei and 62 so that {e^, 61,62} is 
an orthonormal basis for M3, and choose 64 . . . e^-i so that {63 . . . , ed-i} is an orthonormal basis 
for 5"^~^. Now take {rrij} = {61,62,64, . . . ,6rf_i}. Let a, 6 take values 1,2 and let a, [3 take values 
4 . . . , (d — 1) and a, (5 take values 3, . . . , d — 1. By axisymmetry we have that Cabcd vanishes if there 
are an odd number of indices of the form a, (3. We also have 

for some quantities a, C-^, b. Now the WAND condition C^upai^m^l^mj = reduces to 

Cdb = -CoaOb^ (125) 



Note that the special case dSs x 5"* can be written in the form (|118|) (with constant Eq and k = 1). 
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The first equation follows from choosing i,j = d, 6 in the WAND condition and the second by 
choosing = a, (3. 

Now, since we have a product metric, Cabcd is fully determined by the Ricci tensor of M3. 
Hence these conditions give conditions on this Ricci tensor. Using the formulae in |12j (and R^/s = 
{d — A)6aj3), we find that the Ricci tensor of M3 must obey 

i?oo = 2, R,i = \R,,6.^i. (126) 
Similarly, the additional condition for £ to be a multiple WAND reduces to 

i?oa = 0. (127) 

Note that these conditions are invariant under cq — > — cq, which implies that if ^ is a multiple WAND 
then so is n. Hence the spacetime is type D or more special. From now on, we assume that ^ is 
indeed a multiple WAND. Note that we can argue as we did in the second paragraph of this section 
to deduce that there is no loss of generality in assuming that cq is axisymmetric, which we shall 
assume henceforth. 

Using capital letters M, N, . . . to denote coordinate indices in M3 , we can summarize the form 
of the 3d Ricci tensor as 

RmN = (a* + 1)UMUN + (/i - l)gMN, (128) 

where = Cq^. This is the Ricci tensor that would arise from a solution of the 3d Einstein 
equations sourced by a perfect fluid with energy density fj, and pressure p = 1 (with SnG^ = 1). The 
contracted Bianchi identity (or stress tensor conservation) gives 

in + l)u ■ Vu^' = 0, (129) 

and 

(/i + 1)6* + u • V/x = 0, (130) 

where the expansion 6 is defined hy 9 = V • u. The first of these equations implies that either 
fi = —1 or u*^ is tangent to affinely parametrized geodesies in M3. In the former case, we have 
Rmn = —"^ghiNi which implies that M3 is locally isometric to AdS^ with unit radius, which implies 
that M3 X S'^~^ is conformally flat, so the physical spacetime is type O. Assume henceforth that this 
is not the case, so 7^ — 1 and li^^ is geodesic in M3. 
The Einstein equation for the physical metric E'^g is 

M^i = E-^R^ + {d- 2)E-^V^'VyE-^ - E-'^V'^ E'^"'^ d^i , (131) 

where V is the covariant derivative with respect to the unphysical metric g, and indices are raised 
and lowered with this metric. The components tangent to the sphere give 

-l—E-'^V'^E'^-'^ = {d- 4)^-2 _ ^_ (132) 

Using this, and ()128p . the components tangent to M3 give 

E-^ [{fi + l)uMUN + (m - (d - ^))9mn] + {d- 2)VmVnE-^ = 0. (133) 

We now act on this with Vp, antisymmetrize on MP, and use the fact that the Riemann tensor in 
3d is determined by the Ricci tensor, which is given by p28|) . This results in the equation 

= (/i + l)V[pUM]UN -{d- + l)V[pE~^gM]N + E~^{V ypii)uM]UN + E'^{fi + 1)B[mp]Un 
+ E-^ifJ' + '^)Bn[pUm] + -£^"^(V[pm)s'a/]jv - (c^ - 2)(^ + l)gN[MUp]U ■ VE'^ 
- {d-2){fi + l)V[pE-\M]UN (134) 
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where 

BmN = ^NUM- (135) 

Contracting with , this equation reduces to (/i + \)B\^mp] = 0. Since we are assuming that + 1 
is not identically zero, we must have 

B[MN] = 0, (136) 
i.e. du = 0, so u^'^ is hypersurface-orthogonal. Define the projector 

hMN = QMN + UMUN (137) 

and now contract (I134p with Hq to get 

= -{d-3){^, + l)V[pE-^hM]Q + E-'{V[p^l)hM]Q + E-\^l + l)BQ[pUM] 

- {d-2){ij + l)hQiMUp]U-VE-\ (138) 

where 

Bmn = h^jh'^BpQ. (139) 

We can define the expansion and shear of the geodesic congruence tangent to u^^ in terms of the 
trace and traceless parts of Bmn- 

Bmn = -^OhMN + o'mn- (140) 
Contracting (jl38p with gives (/i + 1)-Bqp oc hpq, hence the congruence is shear-free: 

<TMN = 0. (141) 

Equation ()138p now reduces to 

X[phM]Q = 0, (142) 

where 

Xp = -{d - 3)(/i + l)VpE~'^ + E-^Vp^j, - ^9E~\n + l)up -{d-2){fj. + l)upu ■ VE'K (143) 

However, contracting (jl42p with reveals that Xp = 0. Decomposing this into a part orthogonal 
to u^^ and a part parallel to u^^ gives 



h'l,VN[E''-'{fi + l))=0, (144) 

e = -2E-^u ■ VE (145) 

where ()130p was used to simplify the second equation. 

Let So be a surface orthogonal to u^^ (recall that is hypersurface orthogonal), let be 
coordinates on Sq. Assign coordinates (T, x*) to the point proper time T along the geodesic tangent 
to u^'^ starting at the point on Sq with coordinates x*. In this chart, the metric is 

ds"^ = -dT^ + hij (T, x)dx'dx^ , (146) 

and u = d/dT. Prom the definition of Bmn and using the fact that the rotation and shear of the 
geodesies vanish, and equation (jl45p we deduce that 

h,,{T,x) = E-^Hij{x), (147) 
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for some 2- metric Hij independent of T. Eliminating 9 between equations ()130p and ()145p gives {/j, + 
1) = f{x)E'^ for some (non-zero) function /. Substituting this into equation (jl44p and integrating 
gives 

for some function g{T). Now contracting ()133|) with u^^u^ gives 

af.^-^ + £^-1 = 0. (149) 
Using ()148p and the freedom to shift T by a constant we can solve to obtain 

^-1 = r(x)"^sinT, (150) 

for some non-zero function r{x). 

Putting everything together, the physical metric is 



sin^r 



ds^ = r{xf ( ) + Hij{x)dx'dxK (151) 



The metric in brackets is the metric of (d— 2)-dimensional de Sitter space. The full metric is invariant 
under the de Sitter isometry group. Hence if we Wick rotate to Euclidean signature then obtain 
a spherically symmetric spacetime so we can apply Birkhoff's theorem to deduce that the above 
metric must be either the Kaluza-Klein bubble spacetime (1122p . or (if r{x) is constant and A > 0) 
dSd-2 X S\ 
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A Curvature tensors for axisymmetric metrics with axisymmetric 
geodesic WAND 



In this Appendix, we record the non-zero components of the Riemann tensor of the metric ([3] 
describing an axisymmetric spacetime with an axisymmetric geodesic WAND, using the null basis 
([28]) (with B = 0) 

Roioi = l(^^ + 2d'^Uy (152) 
i?0202 = (154) 

RoaOP = -Saf}-^, (155) 

i?ioi2 = [2CDd'^U - DUd'^C - 2Ddrd,U - 2Dd„drC - 2{CdrU - d^U)drD 

+Ad^CdrD - 2drCdyD + UdrCdrD] , (156) 
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-Ri2i2 = ^ [2D(C'^dlU + dlU) - 2CDUd^C - D^{4d^D + U'^d^D) - 4CDdrd,U - AD'^UdydrD 

- 2D{2CdydrC - 2d^d^C - Udrd^C) - 2C^drUdrD + 2CUdrCdrD - 2d^Ud^D 

- 2D^{d^UdrD - drUdyD) - 2D{drUd^C - d^UdrC) + 2C{drUd^D + d^UdrD) 
+AdyC{CdrD - d^D) - 2UdrCd^D] , (157) 



Riaip = -^^^ [-D'^{4.dlE + U^d^E + AUdydrE) + 2C'^drUdrE - 2CUdrCdrE + 2d^Ud^E 

- 2D^{dyUdrE - drUdyE) + 4.d^C{d^E - CdrE) + 2{UdrC - CdrU)d^E 
-2Cd^UdrE] , (158) 



-R2021 = ^ [-2D{Cd^C + DUdlD) - AD^d^drD + 2Ddrd,C 

-D{drCf + 2drC{CdrD - d^D) - 2D'^drUdrD] , (159) 



R2a2(3 = [-D{C^d^E + dlE) + 2CDdrd,E - CDdrCdrE + (C^ - D^U)drDdrE 

+d^Dd^E + {Dd^C - Cd^D)drE - D'^{d^DdrE + drDdyE) - CdrDd^E] , (160) 

Raopi = ^ [D^iUd^E + 2dydrE + drUdrE) + CdrCdrE - drCd:,E] , (161) 

Rocm = [2D{d^EC - drd^E) + {DdrC - 2CdrD)drE + 2drDd,E] , (162) 

Ram = 77^ [2CDUd^E + ACDdydrE - ADdyd^E - 2DUdrd^E - 2DdrCdyE + {2CDdrU 

4i/ ±!j 

-2Dd^U - DUdrC - ACdyD - 2CUdrD)drE + 2{2dyD + UdrD)d^E] , (163) 



Rap^S = [D^ - {C' + D^U){drEf - {d,Ef - 2drE{D^dyE - Cd.E)] . (164) 

The non-zero components of the Ricci tensor are 

i?oo = -^-(d-3)^, (165) 

i?oi = --1^ [D^dlU + CDdlC + D'^Ud^D + 2D'^dydrD - Ddrd^C 
+D{drCf - drC{CdrD - d^D) + D'^drUdrD] 
- ^ [D^Ud^E + 2D^dydrE + drC{CdrE - d^E) + D^drUdrE] , (166) 
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02 



1 

2D 



dr-D 



+2{d-3) 



D 

drDd^E 
DE 



{d-3) 



drE 



E 
drdzE 
E 



drC + 2{d - 3) 



d^E drDdrE 



E 



DE 



C 



(167) 



= [2D{C'^dlU + dlU) - 2CDUdf.C - D^{4d^D + U'^d'^D) - ACDdrd.U - AD^Ud^drD 

- 2D{2CdydrC - 2d^d,C - UdrAC) - 2C^drUdrD + 2CUdrCdrD - 2dzUd,D 

- 2D^{d,jUdrD - drUdyD) - 2D{drUdzC - d.UdrC) + 2C{drUd,D + d.UdrD) 
+Ad^C{CdrD - d,D) - 2UdrCd,D] 

+ ~ \-D^{4d^E + U^d^E + AUd^drE) + 2C'^drUdrE - 2CUdrCdrE + 2d^Ud^E 
4D^E '- 

- 2D^{dyUdrE - drUdyE) + 4d^C{dzE - CdrE) + 2iUdrC - CdrU)d^E 
-2Cd,UdrE] , (168) 



i?i2 = ^ [2CDd^U - DUdlC - 2Ddrd,U - 2Dd^drC - 2{CdrU - d,U)drD 
+Adi,CdrD - 2drCd^D + UdrCdrD] 

+ ^'!~^} \2CDUd'^E + ACDdydrE - ADd^d.E - 2DUdrd,E - 2DdrCd^E 
AD^E '- 

+ {2CDdrU - 2Dd,U - DUdrC - ACd^D - 2CUdrD)drE + 2{2d^D + UdrD)d,E] , 

(169) 

R22 = \-2D{CdlC + DUd^,D) - AD^d^drD + 2Ddrd,C 
2D'^ '- 

-D{drCf + 2drC{CdrD - d^D) - 2D^drUdrD] 

+ ^Ir^ \-D(C^d^E + d^E) + 2CDdrd,E - CDdrCdrE + (C^ - D^U)drDdrE 
D-^E 

+d,Dd,E + {Dd,C - Cd,D)drE - D^{dyDdrE + drDd^E) - CdrDd.E] , (170) 



Ro,p = 5ap [-D{C^ + D'^U)dlE - DdlE - 2D^d^drE + 2CDdrd,E - D^drUdrE 

- 2CDdrCdrE + (C^ - D'^U)drDdrE + d^Dd^E - D'^{d^DdrE + drDd^E) 
+D{drCd,E + dzCdrE) - C{drDd,E + d^DdrE)] 

+ ^^2^ - + D'U){drEf - {d^Ef - 2drE{D^d,E - Cd-^E)]^ . (171) 
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